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© Method and apparatus for measurement of fluid flow in a drilling rig return line. 



© An improved method and apparatus is disclosed 
for use in a return line system (26) of a drilling rig to 
determine a volumetric fluid flow rate and, therefore, 
to aid in detecting a kick in a well. The apparatus 
comprises an ultrasonic level sensor (10), an ultra- 
sonic doppler velocity sensor (12) and a digital com- 
puter (14) electrically coupled to the level sensor 
and the doppler velocity sensor to determine the 
fluid flow regime in the return line system (26). The 
^ flow regime is determined by the programmed com- 
^puter (14) to accurately determine the volumetric 
fluid flow rate in the return line system in real time, 
during drilling operations. 
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METHOD AND APPARATUS FOR MEASUREMENT OF FLUID FLOW IN A DRILLING RIG RETURN LINE 



BACKGROUND OF THE INVENTION 

Reld of the Invention 

This invention relates generally to the fluid of 
measuring the volumetric flow rate of a fluid. More 
particularly, the invention relates to a method and 
apparatus for measuring the volumetric flow rate of 
a fluid in a drilling rig return line in real time. 



FLOW RATE CHANGE 

- An increase in the flow-out or flow rate leav- 
ing the well while pumping at a constant rate is one 

5 of the primary kick indicators. The increased flow 
rate is interpreted to mean that the formation is 
forcing formation fluids into the well bore. A de- 
crease in the flow rate exiting from the well while 
pumping at a constant rate is an indicator of lost 

70 circulation. 



Description of the Prior Art 

Well control and blowout prevention has be- 
come an important concern in the oil and gas 
drilling industry for a number of reasons. Well 
blowouts cause higher drilling costs, possible loss 
of life, and waste of natural resources. An addi- 
tional reason for concern is the increasing number 
of governmental regulations and restrictions being 
placed on offshore drilling operations partially as a 
result of recent, much-publicized well control in- 
cidents. 

A kick can be defined as a well control prob- 
lem in which the pressure found within the drilled 
formation is greater than the mud or fluid hydro- 
static pressure acting on the borehole or face of 
the formation. This formation pressure causes 
fluids to flow from the formation into the well bore. 
In almost all drilling operations, the operator at- 
tempts to maintain a hydrostatic pressure greater 
than the formation pressure and thus prevent kicks. 
On occasion however, and for various reasons, the 
formation pressure exceeds the mud pressure and 
a kick will occur. Kicks have become even more 
common due to the present trend of increasing 
drilling rates by using lighter drilling mud. 

Another problem encountered when drilling a 
well is drilling fluid loss into the formation. This 
problem, known by the shorthand term, "Lost Cir- 
culation", occurs where the drilling fluid is flowing 
into a subterranean formation through which the 
borehole passes. Such condition should be de- 
tected quickly by a driller to prevent damage to 
such a formation and excessive loss of the drilling 
fluid. 

A number of kick or lost circulation "indicators" 
can be observed at the surface before a kick has 
had time to result in a dangerous blowout or exces- 
sive time has elapsed since the beginning of lost 
circulation/Three of these are: 



FLOWING WELL WITH PUMPS OFF 

15 - When the rig pumps are not moving the mud, 
a continued flow-out from the well indicates that a 
kick is in progress. An exception to this indicator is 
when the mud in the drill pipe is considerably 
heavier than that in the annulus, as in the case of a 

20 slug. 



PIT VOLUME CHANGE 

25 - If the volume of fluid in the pits is not 
changed as a result of surface controlled actions, 
an increase in pit volume indicates that a kick is 
occurring. The fluids entering the well bore as a 
result of the kick displace an equal volume of mud 

30 at the flow line and result in a pit gain. A decrease 
in pit volume under these conditions indicates lost 
circulation. 

Two of the kick early warning signs described 
above require measurement of an increase in flow 

35 rate from the fluid return line while the other re- 
quires measurement of an increase in pit volume. 
These indicators are difficult to interpret when drill- 
ing from a floating drilling vessel because of the 
heaving and rolling of the drilling vessel in re- 

40 sponse to wind and waves. Floating drilling vessel 
heaving and rolling creates fluid return line flow 
rate changes. 

It has been found that the time elapsed be- 
tween the beginning of a kick.deep in the well and 

45 its detection at the surface by pit level monitoring 
is too long to provide sufficient time to bring the 
well under control, usually by adding weight to the 
drilling fluid. 

Studies have shown that accurate differential 

50 flow measurement, of the order of 25 gallons per 
minute (GPM) of a total of 1200 GPM, provides the 
earliest possible surface detection of kicks and/or 
lost circulation. This method requires the use of 
flow meters of high absolute acuracy under widely 
varying conditions, for both flow-in and flow-out 



systems. 

Presently, flow-in measurement is based on the 
number of strokes per minute of triplex mud pumps 
(see Figure 1). The flow rate obtained from the 
pump strokes is then corrected by a volumetric 
pump efficiency. This pump efficiency can fluctuate 
between 80% to 95%, accounting for inaccuracies 
of ± 10% in the flow-in measurement. 

The prior flow-out measurement has usually 
included a "paddle" system installed in the rig 
return line. The paddle is a hybrid flow meter 
based on level and target (force) measurements. 
The prior art paddle has an uncaJibrated accuracy 
of around 40%. With calibration on the rig site, the 
"absolute" flow-out measurement is still only ac- 
curate to 10% or 15%, due to the basic non- 
linearity of the device, and due to very poor zero 
stability of the device. Poor zero stability requires 
frequent re-calibration. 

Using the pump strokes and the paddle mea- 
surement for flow-in and flow-out, respectively, the 
best accuracy for the differential flow over the 
entire fluid flow range, cannot be better than 25% 
or 300 GPM (gallons per minute) in 1200 GPM. 
This is more than 10 times the required accuracy, 
rendering prior methods of differential flow rate 
measurement inadequate for desired kick detec- 
tion. 

Electro-magnetic flow meters have also been 
used but have drawbacks. They do not work in oil 
based muds (conductivity too low). They require 
complete modification of the return line. In offshore 
rigs where the return line is short and space is 
limited, there is usually no way to install them. 
They require expensive maintenance to sustain 
their accuracy. 

Other industries have developed flow measur- 
ing systems suitable for use in waste water moni- 
toring systems where sewer outflows must be mon- 
itored for pollution control purposes. These sys- 
tems obtain flow measurements based on the ve- 
locity of the fluid in a channel and the area of the 
channel occupied by the flowing liquid. Ultrasonic 
level detectors and doppler type velocity detection 
units have been used for these applications. 

U.S. Patent 4,217,777 to Newman issued Au- 
gust 19. 1980 discloses such a system and is 
incorporated herein for essential material and for all 
other purposes. Also U.S. Patent 4,202,211 to Per- 
ry issued May 13, 1980 discloses a similar system 
and is incorporated herein for essential material 
and for all other purposes. 

Ultrasonic level detection systems are known in 
the art. Such systems are described in U.S. Patent 
No. 4,024,766 to Perry issued May 24, 1977 and in 
U.S. Patent No. 4,145,914 to Newman issued 
March 27. 1979, both of which are incorporated 
herein for ail purposes. 



The fluid flow in the rig return line is character- 
ized as supercritical flow. None of the known prior 
art accurately measures volumetric flow rate for 
flows in supercritical flow conditions. Though the 

5 above four U.S. Patents disclose ultrasonic measur- 
ing devices, they do not disclose or teach distin- 
guishing between flow regimes which include 
laminar and turbulent flow so as to determine suf- 
ficently accurately a volumetric fluid flow rate which 

to can be used for real time kick or lost circulation 
detection on a drilling rig. 



Identifications of Objects of the Invention 

76 

It is an object of the invention to provide an 
improved method and apparatus for the measure- 
ment of volumetric fluid flow rate which Is signifi- 
cantly more accurate than prior art methods and 
20 provides the measurement in real time while drill- 
ing. 

It is another object of the invention to provide 
an improved fluid flow rate measuring system 
which obtains volumetric flow measurements based 

26 on the measured velocity of the fluid in a conduit 
the cross-sectional area of the conduit occupied by 
the flowing fluid, and a determination as to whether 
the flow is laminar or turbulent. 

It is a still further object of the invention to 

30 provide an improved system for accurately measur- 
ing flow rate by distinguishing between flow re- 
gimes. 

It is still another object of the invention to 
provide an improved method and apparatus for use 
36 in a return line of a drilling rig to quickly and 
accurately detect a kick or lost circulation in the 
well bore. 

It is another object of the invention to provide a 
volumetric flow measuring system for accurately 
40 measuring flow rate of a fluid in supercritical flow 
conditions. 

It is still another object of the invention to 
provide all the foregoing measurements on a float- 
ing rig while simultaneously suitably compensating 
46 for the effects of rig heave. 



SUMMARY OF THE INVENTION 

50 According to the invention, an improved meth- 
od and apparatus is disclosed for use in a return 
line system of a drilling rig to aid in accurately 
determining a volumetric fluid flow rate and, there- 
from, to quickly and accurately detect a kick or lost 

55 circulation in the well bore. 

The apparatus includes an ultrasonic level sen- 
sor, an ultrasonic doppler velocfty probe and a 
digital computer electrically coupled to the level 
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sensor and the doppier probe to determine, under 
program control, the fluid flow regime in the return 
line system. The flow regime is used by the com- 
puter to accurately determine the volumetric fluid 
flow rate in the return line in real time. 

The level measurement is made with an ultra- 
sonic level sensor through the air above the flowing 
fluid. The velocity is measured by an ultrasonic 
doppier probe which takes its measurement in a 
well defined volume of fluid at the location of the 
level measurement. 

The digital computer under program control 
processes input data in conjunction with stored 
input data to determine whether the flow in the 
return line is laminar or turbulent In order to obtain 
an accurate determination of level of the flowing 
fluid, the output of the level sensor is temperature 
compensated by means of temperature data mea- 
sured at two, preferably three locations representa- 
tive of the temperature of the air through which 
ultrasonic pulses travel. Compensation is neces- 
sary because the mud is hot (180°F) and the 
stop/start nature of drilling does not allow thermal 
equilibrium to be established. 

In order to increase accuracy of the flow rate 
measurement, any sediment level in the return flow 
line is measured under zero flow conditions with 
the level sensor. An accurate determination of sedi- 
ment level insures that an accurate determination 
of the cross-sectional area of flow is determined 
which is used in the flow rate calculation. 

The flow regime is determined by first assum- 
ing that the flow is laminar. That assumption is 
tested periodically with real time measurement of 
level h and velocity v with stored numerical in- 
formation of the geometry of the return line. A 
correction factor for laminar flow Is determined 
which when multiplied by the measured velocity 
and the measured cross-sectional area of the re- 
turn line (as a function of measured level) yields an 
accurate determination of instantaneous laminar 
flow rate. 

A value of viscosity associated with such 
laminar flow is determined and then the Reynolds 
number associated with such viscosity is calcu- 
lated. A test is then made to determine if such 
Reynolds number determination satisfies the 
laminar flow assumption. 

If the laminar flow assumption is not satisfied 
by the Reynolds number test, the flow regime must 
be turbulent. By means of instantaneous velocity 
and level measurement information, the local fric- 
tion factor at the measurement location is deter- 
mined using an analysis comparing free gravity 
flow to actual measurements. The local friction 
factor is then used to obtain an accurate correction 
factor for the turbulent flow velocity profile. Such 
correction factor is multiplied with the instanta- 



neous velocity and the instantaneous cross-sec- 
tional area to yield an accurate determination of 
turbulent regime flow rate. 

An averaging process which uses a variable 

s averaging period is used to determine an average 
flow rate for floating rigs where the instantaneous 
flow rates are varying as a result of the heave and 
pitch of the rig. 

A method of using the level sensor alone is 

io presented based on stored data when both the 
velocity and level sensors were known to be op- 
erating properly. 

Output measurements determined periodically 
of average output flow rate, differential flow rate 

is between the output flow from the well and the flow 
into the well, and viscosity as determined from the 
flow regime calculations are recorded as a function 
of time. The differential flow rate is used to drive a 
driller's gauge as a warning of a kick or lost circula- 

20 tion in the well. 



BRIEF DESCRIPTION OF THE DRAWINGS 

25 The objects, advantages and features of the 
invention will become more apparent by reference 
to the drawings which are appended hereto and 
wherein like numerals indicate like parts and 
wherein an illustrative embodiment of the invention 

30 is shown, of which: 

Figure 1 is a schematic view of a drilling 
fluid or "mud" circulation system for a floating 
drilling rig where a flow measuring system em- 
bodying the invention may be used; 

35 Figure 2 is a cross-sectional view of a return 

line system of the invention illustrating the retrofit- 
ting of an ultrasonic sensor system into a preexist- 
ing opening for a prior art paddle in the return line 
and schematically illustrating connection of the 

40 sensor system to a computer and associated in- 
strumentation; 

Figure 3 is a section view along line 3-3 of 
Figure 2 illustrating the level sensor location to 
determine the cross-sectional area of fluid; 

45 Figures 4A and 4B present a logic flow dia- 

gram illustrating the method of the system; 

Figure 4C illustrates a laminar flow pattern 
and shows the velocity sensor measure point in the 
fluid and is used to illustrate determination of the 

so laminar conversion constant used to determine 
average velocity for laminar flow; 

Figure 4D illustrates a turbulent flow pattern 
and shows the use of the f^ factor in determining 
the initial slope of the turbulent flow profile so as to 

55 more accurately determine the turbulent conversion 
constant used to determine average velocity for 
turbulent flow; 



4 



Figures 5A and 5B illustrate respectively a 
cross-section of a return line with temperature sen- 
sors located at top, bottom and mid-way positions 
between a level sensor and the fluid and a graphi- 
cal example of the determination of a functional 
* relationship between temperature and the distance 

between the head of the level sensor and the fluid 
in the return pipe; 

Figure 6A illustrates the relationship between 
laminar flow rate and a parameter h/D used in the 
determination of viscosity and whether the flow is 
laminar or turbulent; 

Figure 6B shows relationships between fric- 
tion factor of the return line as a function of Reyn- 
olds number used in the determination of whether 
the flow is laminar or turbulent; 

Figure 6C shows the relationship between 
the parameter h/D and Cv l for laminar flow, where 
Cvl is the laminar velocity profile correction factor 
used to determine average velocity for laminar 
flow; 

Figure 7 shows a relationship between the 
average friction factor f^ for the return line and 
the local friction factor f^ at the position in the 
return line where the level and velocity are mea- 
sured; and 

Figure 8 illustrates the method by which 
average flow rate is determined from instantaneous 
flow rates under conditions of a heaving drilling rig. 



DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

The volumetric flow sensor system, generally 
designated S, and associated computer 14 em- 
bodying the invention is illustrated in Figures 1-4. 
The flow measuring system includes an ultrasonic 
level sensor 10, an ultrasonic doppler velocity 
probe 12 and a processing means or digital com- 
puter 14, as best illustrated in Figures 2 and 4. 

Referring to Figure 1 , a floating drilling rig mud 
circulation system is schematically illustrated. Al- 
though a floating drilling rig is illustrated, the inven- 
tion may also be used with a bottom supported 
offshore drilling rig or a land drilling rig. As dis- 
cussed above, flow rate into the well may be de- 
rived by counting the strokes per minute of mud 
pumps 16. After the "mud" or drilling fluid travels 
down the drill string 18, it moves up the annulus 20 
between the casing 22 and the drill string 18 to bell 
nipple 24. A return line 26 communicated with the 
bell nipple 24, as best shown in Figures 1 and 2, to 
return the mud to the mud pits 28. The flow-out 
measuring sensor system S according to the inven- 
tion is disposed in the return line 26 at a lower 
elevation than the bell nipple 24. 

The sensor system S in cooperation with com- 



puter 14 is provided to accurately measure In real 
time the flow Q through return line 26. Volumetric 
flow Q is simply the product of the cross-sectional 
area A of fluid flowing at a given location in the line 
s and the average velocity of the fluid moving at that 
location. 

The sensor system S and computer 14 ac- 
curately determine the cross-sectional area A of 
pipe 26 through measurement of the level L from 

to the top of the line 28 to the surface of the drilling 
fluid. With a measurement of the sediment level in 
the line 26 when no fluid is flowing in the line, an 
accurate determination of the area A of the flow 
may be determined In real time. 

75 Measurement of the velocity at a point in the 
fluid flow stream, with knowledge of the geometry 
of the return line allows a determination of the flow 
regime, whether turbulent or laminar, and conse- 
quently allows a more accurate determination of 

20 average velocity of the flow. Accurate determina- 
tion of cross-sectional area and average velocity 
results in a more accurate determination of drilling 
return line flow rate, and when combined with flow- 
in measurements from the mud pumps, results in 

25 an accurate timely indicator of a "kick" or lost 
circulation in the well. 

The return line 26 usually has a circular cross- 
sectional area, as best shown in Figure 3. Other 
cross-sectional profiles, for example, elliptical, rec- 

30 tangular or square, could be used. For a preferred 
circular return line 26, the internal diameter D of 
the return line 26 and the fluid height h allows 
determination of the fluid cross-sectional area A. 
Figure 3 indicates that ft is necessary to determine 

35 if a sediment level exists in order to make an 
accurate determination of the area A of the flowing 
fluid. Measurement of sediment level SED (if any) 
is described in more detail below. 

As illustrated In Figure 2, sensor system S is 

40 preferably constructed to fit upon the flange 11 of 
an opening in drilling line 28 which has been pro- 
vided for prior art paddle velocity sensors. It of 
course may be provided in any other top entry 
opening in drilling line 26. Thus sensor system S 

45 achieves access to the interior of line 26 without 
the need for line 26 being disconnected from an 
existing rig hook-up to bell nipple 24. 

Alternatively, probe 52 may be inserted 
through a hole in the bottom of Une 26 in the 

so vicinity of the level sensor 10 so that there is no 
obstruction in the line 26. Such a placement of 
probe 52 may be indicated where "gumbo" clay 
formations are expected to be encountered during 
drilling. 

55 The sensor system S includes a level sensor 
10 for measuring the distance L + L', and from that 
measurement, information for the determination of 
h (see Figure 3). Level sensor 10 includes a hous- 
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ing 40 in which an ultrasonic transmitter/receiver is 
mounted. Ultrasonic pulses from a transmitter with- 
in sensor 42 are directed toward the fluid F and 
return to a detector within sensor 42 to produce a 
signal on lead 44 proportional to L + L'. 

The preferred level sensor 10 includes an ultra- 
sonic head 42 and a microprocessor 45. The head 
42 includes a single crystal transmitter/detector and 
is constructed of PVC with a Buna-N acoustic win- 
dow. A sound cone dispersion of seven degrees is 
produced by the crystal. The crystal is a low Im- 
pedance, passive capacitative element driven by a 
driving voltage of about 350V (RMS) at about 40 
KHz at a 10% duty cycle. 

The microprocessor 45 controls the head 42 
and receives differential time information between a 
generated pulse and its echo from the fluid. An 
analog signal on lead 44 proportional to the dis- 
tance L + L' from the head 42 to the fluid (see 
Figure 3) is applied to computer 14 via conven- 
tional A/D converters 47 for further processing. A 
temperature sensor 46, disposed near the head 42 
applies an analog signal via lead 48 to A/D convert- 
ers 47 for further processing discussed below. 
Computer 14 is preferably a Digital PDP-11 
minicomputer. 

A temperature sensor 53 associated with ve- 
locity probe 52 produces a signal representative of 
temperature at the fluid level on lead 43 to micro- 
processor 41. An analog signal representative of 
fluid temperature on lead 43' is applied from micro- 
processor 41 to A/D converters 47 for application to 
computer 14, 

An additional temperature sensor 29 may ad- 
vantageously be placed on flange 13 such that a 
temperature measurement is made at the disconti- 
nuity between housing 40 and return line 26 is 
applied to computer 14 via lead 27 and A/D con- 
verters 47. 

In a preferred embodiment of housing 40, an 
extension 50 of the housing 40 extends partially 
into the interior of fine 26 so as to prevent exces- 
sive angular dispersion of the pulses and false 
echo returns from the discontinuity. 

The velocity probe sensor 12 is preferably an 
ultrasonic doppler velocity sensor known in the art 
of sonic velocity sensors. Typically, such sensors 
contain two piezo-eiectric crystals mounted at the 
end of the probe 52 in the vicinity of the height 
measurement by head 42. One crystal continuously 
emits ultrasonic waves into the fluid, and the other 
crystal receives the ultrasonic waves reflected by 
the fluid particles. These received waves have a 
shift in frequency with respect to the transmitted 
waves proportional to the velocity of the reflecting 
particles. A doppler shift signal representative of 
the velocity measured by probe 52 in the fluid is 
applied to microprocessor 41 via lead 39. Micro- 



processor 41 preferably performs digital processing 
of the doppler signals. An analog output signal 
proportional to the velocity of the fluid at probe end 
52 is applied on lead 39' to A/D converters 47 for 

5 input to computer 14. 

The measuring tip or end 52 of velocity sensor 
12 is preferably disposed within the return line 26 
to give a fixed measuring point approximately one- 
third the maximum fluid flow height from the bot- 

10 torn of the line 26. This fixed measuring point is 
entered into the computer 14 for calculation of the 
average velocity from the measured velocity for 
both turbulent flow and laminar flow. 

A wash head 56 is advantageously provided 

75 within sensor S to wash probe end 52 of sensor 12 
during zero flow conditions. Pressurized water 
source 58 provides the water to accomplish clean- 
ing of the probe thereby obviating the need for 
manual cleaning of the probe if mudcake buiid-up 

20 were to occur. Wash head 56 is also used in the 
zero flow measurement described below to deter- 
mine whether or not velocity probe sensor 12 is 
properly working. 

25 

Computer Processing Of Input Information 

The processing of level and velocity informa- 
tion to produce an accurate determination of return 

30 line volumetric flow rate may be better understood 
from the flow chart of Figures 4A and 4B. The flow 
chart is a logic flow diagram having logic boxes 
representative of a digital program stored in ROM 
memory 34 (Figure 2) of computer 14. The com- 

35 puter 14 is reset at the beginning of the operation 
which automatically establishes registers in the 
CPU 32 to hold the address of memory locations in 
RAM 30 that are to be accessed for data. Such 
data include the stored input parameters as in- 

40 dicated in box 100. 

As indicated by box 200, signal input measure- 
ments via leads 64, 44, 48, 27, 39, 43, 92 and 96 
are continuously monitored by computer 14. CPU 
32 under program control periodically performs the 

45 steps of Figures 4A and 4B to produce output 
signals representative of Q out , the flow rate of pipe 
26, AQ the difference between Qo Ul and Q|„, the 
input flow rate into the well, and jjl, the viscosity of 
the flowing fluid. 

so 

Temperature compensation of AT measurements 
to determine instantaneous height of flow 

55 A determination of height h is periodically 
made as indicated in box 106 with input information 
of transit time AT from level detector sensor 10 
and the temperature measurements along the path 
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of the ultrasonic pulse between the head 42 and 
the fluid flowing in return pipe 26. The distance M 
» L + L' to the top of the returning drilling fluid is 
computed in calculating unit 28 in terms of the 
product of the velocity of sound and transit time 
AT. The velocity of sound is preferably compen- 
sated by the temperature measurements, tu, U and 
U. 

Figures 5A and 5B illustrate the technique for 
compensation. A relationship t » f(x) between ab- 
solute temperature and the distance between the 
fluid and the ultrasonic head is determined, for 
example by fitting a polynomial between the three 
points x = 0 f t=tu; x=x„ t=U; and x = x„ t-tl. Of 
course if only the temperature at the head, t„ and 
the fluid U are available, a straight line may be 
determined as the functional relationship between t 
and x. 

The velocity of sound in the air through which 
the ultrasonic pulses are passing is determined by 
calculating an average velocity of sound: 



Determination of sediment level (SEP) 

The level sensor 10 is preferably instructed 
under program control by the computer 14 to pro- 
vide data at zero flow for sediment in the return 
line. Conventional flow indicating means in the re- 
turn line may alternatively be electrically coupled to 
the computer 14 to provide the data of zero flow. 

The apparatus illustrated in Figure 2 may be 
used to measure sediment level in the return line at 
zero flow, but a determination as to the proper 
operation of velocity probe 52 is desirable. Box 103 



of Figure 4A illustrates this periodic determination 
under computer control. To achieve this objective, 
a spray jet of water from head 56 from water 
source 58 is directed under program control to 
s probe 52. If a velocity is measured from probe 52, 
via microprocessor 41, then logic stored in pro- 
gram control of computer 14 determines that the 
velocity information from probe 52 is correct and 
that the level information from head 10 and micro- 
10 processor 45 is correctly measuring sediment level 
SED. Logic box 105 illustrates the determination of 
sediment level (SED). 

The distance to the top of the sediment is 
computed in the calculating unit 28 in terms of the 
is product of the velocity of sound and transit time 
AT. The level of the drilling fluid (h) is then com- 
puted by subtracting the sediment distance SED 
from the calculation described above, that is: h * 
D-<M-L>SED. The determination of sediment 
20 and/or foreign matter In the return line provides a 
more accurate cross-sectional area A as shown by 
conventionai determination of A with reference to 
the geometry illustrated in Figure 3 and, therefore, 
a more accurate volumetric flow rate is determined. 



Determination of flow regime 

As indicated by box 108 in Figure 4B, a deter- 
mination of the flow regime, that Is whether the flow 
is laminar or turbulent is made under program 
control The determination begins with an assump- 
tion that the flow is laminar. This assumption is 
tested periodically in real time with input data of 
velocity v and level h measured by sensor S. It is 
known that the flow rate for laminar flow is a 
function of the kinematic viscosity of the fluid 

U 
0 

the level of the flowing fluid h, its velocity v and the 
return line geometry, that is, 
Qlam = f<*. h, v, D, e, g). 
As illustrated by Figures 4C and 6C, an exact 
correction factor Cv l can be computed for laminar 
flow because of Its parabolic velocity profile and 
from information of the measured level h and the 
height or position P where velocity v is measured. 
Thus, the flow, Q tem is, 

Qlam = VX AX Cvl (1) 

based on the area determination of box 104, as 
shown in Figure 4A and the velocity signal mea- 
surement. 

Next under program control, a tabular relation- 
ship between 



1 X"X 

C - — / 2 / YRf (x) dx 2s 
avg X 2 x=0 25 



30 



where y and R are predetermined constants and f- 
(x) is the functional relation obtained by fitting a 
curve between the absolute temperatures mea- 
sured at the distances between the fluid and the 
head. Such measurements and calculation insure 
an accurate determination of level h and area A of 
the fluid flowing in the line. The distance M, equal ^ 
to L+L' as illustrated in Figure 3 is then deter- 
mined as 

M = L + L' = AT(C 8VQ ). 
The diameter D of the cylindrical return line 26 
and the distance L' from the head 42 to the return 
line 26 are obtained from the stored parameters in 
memory 34. The level of the fluid is then deter- 
mined: 

h = D-(M-L). 
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2 Q lam v 



(g sin 8) D 4 



and h/D 



is addressed in computer 14 as illustrated in Figure 
6A. That is, with knowledge of the measured level 
h, and hence h/D, the corresponding tabulated val- 
ue of 



2 Q 



lam v 



gfsin e)o 4 



is determined, and hence a value of ?* is deter- 
mined knowing Q| 8mf (g sin B) and D 4 , as defined 
by box 100. 

Next, under program control, the value of / is 
used to calculate a Reynolds number R for the 
flow: 



where, as before, v is the measured velocity, and 
d h is the hydraulic diameter of the pipe, 

A 

,a h -— > 

where A is the cross-sectional area of the fluid flow 
as determined in box 104 and P is the "wetted" 
perimeter of the actual fluid flowing (calculated 
from D and h). 

As illustrated in Figure 6B, if R* as determined 
from the calculations above is less than 2100, the 
flow is actually laminar and the flow Q tam calculated 
above (Equation 1) is the correct value of the flow 
rate in the pipe 26. If that determination is made by 
the computer program in logic box 108, the correc- 
tion factor Cyi is as previously determined (Logic 
Box 110) from the geometrical relationship between 
C v and h/D as illustrated in Figure 6C. Having 
determined that the flow is laminar, the kinematic 
viscosity v is also correct and the viscosity u = v p 
can be found from the density p. 

If R* is greater than 2100, the flow regime must 
be turbulent flow as illustrated in Figure 6B. With 
turbulent flow, the velocity profile is almost uniform, 
but to determine a very accurate value for the 
correction factor to be applied to the turbulent flow 
velocity, the following sequence of calculations is 
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made in logic box 111. 

First a turbulent flow rate is determined by 
assuming that Cv « 1:Q = Axvx (1). 

Next, the energy Ec of the fluid flowing under 
critical flow conditions is calculated: 




where Q is the flow rate assumed from step 1 . 

Next the energy loss hf in the return line 26 is 
calculated: 



2g 



where Z, and g are as defined in stored parameter 
box 100 and h and v are the measured level and 
velocity. 

Next, the average friction factor over the whole 
lateral flow length LEL = (H* + Z 2 ) 1 * of the pipe is 
determined: 



hf 



"avg 




Finally a predetermined relationship f Ioca , = F- 
(favg) between the friction factor f loca , at the point of 
the level and velocity measurement to f avg , which is 
stored in memory of computer 14, is addressed to 
yield fio C ai. 

Figure 7 graphically illustrates the relationship 
between f toca i and favg- The correction factor Cvt is 
determined more accurately having a value of local 
friction factor, f loca!r in that the friction factor at the 
local measuring point for level h and velocity deter- 
mines the initial slope of the turbulent velocity 
profile. Figure 4D illustrates the relationship be- 
tween velocity and level of flow at the velocity 
measurement point P. The value of fiocai defines the 
initial slope for the velocity profile for turbulent flow 
and allows a more accurate determination of Cyt. 

The calculations described above are carried 
out under program control represented by box 111. 
With the correction factor Cyt, the instantaneous 
flow rate for turbulent flow is calculated as Q ou t = 



8 



A x v x 0*. 

Viscosity of the turbulent flowing fluid is deter- 
mined from the relationship between f avg and the 
Reynolds number for turbulent flow as illustrated in 
Figure 6B. 



Determination of Average Row Rate for a Floating 
Drilling Rig 

The instantaneous flow rate determinations as 
described above may be used directly where the 
apparatus of Figure 2 is used on a bottom sup- 
ported marine drilling rig or a land rig, but where 
the rig is a "floater", the instantaneous flow rates 
vary with the heave and pitch of the vessel. For 
that reason as "average" flow rate must be deter- 
mined for an accurate determination of CW in 
order to detect a dangerous kick or fluid loss 
condition in the well in combination with Q| n data 
from mud pumps 16. 

Figure 8 illustrates a preferred method for aver- 
aging instantaous flow rate data Q output from 
logic box 112. The position or "height" of the rig as 
a function of time is measured and stored as a 
function of time from sensor 90 via lead 92 to 
computer 14 along with vector information as to 
whether the position has increased or decreased 
from a previous time increment At. The instanta- 
neous flow rate Q t is associated with that position 
and vector as a function of time. For example, 
Figure 8 shows that the t = t* an instantaneous 
value of Q t and the rig's position P, and up vector t 
is stored. 

The averaging method determines an average 
Q at each time interval, for example t = t«. The 
method, according to the invention, is to look back- 
ward in time over the stored data to find a previous 
position with the same vector. For example as 
illustrated by Figure 8, the first time in which the 
position and vector is the same as that of time t ■ 
tt is at t = t 2 . Selecting the period t = tt backward 
to t * t, defines the averaging period T which is 
used to average the instantaneous values of flow 
rate Q associated with each time increment. 

The averaging method described above is ad- 
vantageous in that it achieves an average relatively 
quickly in terms of the time period T of the actual 
up-down heave motions of a floating drilling rig, 
and over the period T the net contribution of heave 
to the flow out is zero. 



Use of Level Sensor Alone to Measure Row 

If as discussed above, the computer 14 under 
program control in logic box 103 determines that 
the velocity sensor is no longer working, control is 



passed to logic box 1 16, as shown in Rgure 4B. 

Logic box 116 indicates that while the velocity 
sensor was previously working, data relating Qout(t) 
and level (h) are being applied to box 116 where 

s they are stored In memory. An empirical relation- 
ship between CWt) and level (h) is thus formed in 
memory based on actual real time operating his- 
tory for the rig and sensors. When logic box 116 is 
entered by virtue of a signal that the velocity sen- 

?o sor Is no longer working, the experience table 
between Q{t) and h(t) is entered with an instanta- 
neous value of h(t) to produce an output of in- 
stantaneous flow rate O^t) which is then passed 
to logic box 114 for the average flow determination 

75 as described above. 



Determination of Differential Row Rate 

20 As illustrated by logic box 118. the difference 
between the average output flow rate, Qoot, as 
output from logic box 114 and the input flow rate 
into the well Q,„ as Input from mud pumps 16 is 
calculated and appHed to recorder 60 and gauge 

26 62. Average flow rate CUxt and viscosity u are also 
applied to recorder 60 for producing a graphical 
record of the flow rate and viscosity of the flow in 
the return line. A negative reading of gauge 62 
indicates "lost circulation" danger, while a positive 

so reading indicates kick danger. 

Advantageously the method and apparatus of 
the invention measures flow rate in a drilHng rig 
return line with improved accuracy. It may be dis- 
posed on the standard connector of the return line 

36 obviating extensive rig modifications for Its deploy- 
ment. It works with water based, or oil based muds 
and may be used under any kind of flow conditions 
such as tranquil or shooting, laminar or turbulent 
The apparatus not only measures flow rate, but 

40 also viscosity of the mud. Since the measurement 
is an absolute flow measurement no field calibra- 
tion is required and can be employed on heaving 
offshore drilling rigs. Averaging of the heave period 
allows a time average of return line flow to be 

46 determined. 

Various modifications and alterations in the de- 
scribed methods and apparatus will be apparent to 
those skilled in the art of the foregoing description 
which does not depart from the spirit of the tnven- 

50 tion. For this reason, these changes are desired to 
be included in the appended claims. The appended 
claims recite the only limitations of the present 
invention and the descriptive manner which is em- 
ployed for setting forth the embodiments and is to 

56 be interpreted as illustrative and not limitative. 
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Claims 

1. A system (S) for determining the volumetric 
flow rate of material (F) flowing in a drilling rig 
return line (26) partially filled with said material, 
said system comprising, 

an ultrasonic pulse echo measurement sensor 
(10) for generating a first signal indicative of ma- 
terial level in said return line; 

an ultrasonic doppler velocity measurement 
sensor (12) for generating a second signal indica- 
tive of the velocity of said material; said system 
being characterized by including: 

a processor (14) responsive to said first and 
second signals for determining the volumetric flow 
rate corrected for the effects of a laminar or turbu- 
lent flow regime. 

2. The system as recited in claim 1 further 
being characterized by including a system (28) for 
determining when there is no material flow in said 
return line and for determining the level of sedi- 
ments in said return line when there is a zero flow 
determination, and a memory (30) for retaining the 
last sediment level determination for use in deter- 
mining volumetric flow rate when material flow re- 
sumes. 

3. The system as recited in claim 1 or 2 further 
characterized by including a system (56) for clean- 
ing said ultrasonic doppler probe sensor (52) when 
there is no flow in the return line. 

4. The system as recited in claim 1, 2 or 3 
further characterized by including a comparator 
(118) for comparing the volumetric flow rate in the 
drilling rig return line (26) with an indication of 
volumetric flow of material flowing into a borehole 
associated with said drilling rig. 

5. The system as recited in one of the 
preceeding claims wherein the system is mounted 
on a floating platform further characterized by in- 
cluding a sensor (90) for measuring the vertical 
position (P) of the floating platform and processing 
apparatus (14) responsive thereto for averaging the 
volumetric flow determination over a variable period 
determined by the time required for the floating 
platform to complete one complete heave cycle. 

6. A method for determining a volumetric flow 
rate of a materia! (F) flowing in an inclined drilling 
rig return line (26) at least partially filled with flow- 
ing material, comprising the steps of: 

measuring the level of said material in said 
return line; 

measuring the flow velocity of said material in 
said return line; 

in response to the level measurement and to a 
predetermined return line geometry, determining 
the cross-sectional area of material flowing in said 
return line; 

said method being characterized by determin- 



ing the average velocity of material flowing in said 
return line corrected for the existing laminar or 
turbulent flow regime in response to the velocity 
measurement, and by determining said volumetric 
5 flow rate in response to the corrected average 
velocity and to said cross-sectional area determina- 
tion. 

7. The method as recited in claim 6 wherein 
said step of determining the cross-sectional area of 

w material flowing in said return line is characterized 
by including the steps of determining when there is 
no flow in said return line and of determining the 
level of sediments on the bottom of said return line 
during the periods of no flow. 

15 8. The method as recited in claim 6 or 7 
wherein said step of determining average velocity 
of material flowing in said return line is character- 
ized by including the steps of measuring a local 
material velocity at a single predetermined local 

20 point in said return line and utilizing a laminar flow 
velocity profile model and said cross-sectional area 
to convert the local velocity into an average ve- 
locity when the flow regime of the flow in the return 
line is laminar. 

25 9. The method as recited in claim 6, 7 or 9 
wherein said step of determining average velocity 
of material flowing in said return line is character- 
ized by including the steps of: 

-in response to the measured velocity, 

30 implementing a free gravity flow model to deter- 
mine the energy loss in said return line (26) and to 
determine a friction factor f as a function of deter- 
mined energy loss of the form, 

35 

f = K | energy loss j 
^ V2 ' 

^ where k is a constant of proportionality and V is the 
measured velocity in said return line; 

-in response to said friction factor f, determin- 
ing the average velocity of said material corrected 
for the effects of turbulence when said flow is 
turbulent and 

45 

-in response to said cross-sectional area and 
said corrected average velocity, determining volu- 
metric flow rate corrected for the turbulent flow 
regime. 

50 10. The method as recited in one of the 
preceeding claims 6 to 9 further characterized by 
including the steps of generating a data base of 
measured level and of volumetric flow rate based 
on level and velocity measurements and accessing 

55 said data base with a level measurement to deter- 
mine volumetric flow rate when a velocity measure- 
ment is not available. 



10 



1 1 . The method as recited in one of claims 6 to 

10 wherein the return line is mounted on a floating 
platform subject to heave and is further character- 
ized by the steps of measuring the vertical position 
(?) of the floating platform and in response thereto, 
averaging the volumetric flow determination over a 
variable period determined by the time period re- 
quired for the floating platform to complete one 
complete heave cycle. 

12. The method as recited in one of claims 6 to 

11 further characterized by comparing the volu- 
metric flow rate in the drilling rig return line (26) 
with an indication of volumetric flow of material 
flowing into a borehole associated with said drilling 
rig. 
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